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Abstract 

The influence of an axial through-flow on the spatiotemporal growth behavior of different vortex 
structures in the Taylor-Couette system with radius ratio rj = 0.5 is determined. The Navier 
Stokes equations (NSE) hnearized around the basic Couette-Poiseuihe flow are solved numericahy 
with a shooting method in a wide range of through-flow strengths Re and different rates of co- 
and counterrotating cylinders for toroidally closed vortices with azimuthal wave number m = 
and for spiral vortex flow with m = ±1. For each of these three different vortex varieties we 
have investigated (i) axially extended vortex structures, (ii) axially localized vortex pulses, and 
(in) vortex fronts. The complex dispersion relations of the linearized NSE for vortex modes with 
the three different m are evaluated for real axial wave numbers for (i) and over the plane of 
complex axial wave numbers for {ii, Hi). We have also determined the Ginzburg-Landau amplitude 
equation (GLE) approximation in order to analyze its predictions for the vortex stuctures {ii, 
Hi). Critical bifurcation thresholds for extended vortex structures are evaluated. The boundaries 
between absolute and convective instability of the basic state for vortex pulses are determined with 
a saddle-point analysis of the dispersion relations. Fit parameters for power-law expansions of 
the boundaries up to Re'^ are listed in two tables. Finally, the linearly selected front behavior of 
growing vortex structures is investigated using saddle-point analyses of the dispersion relations of 
NSE and GLE. For the two front intensity profiles (increasing in positive or negative axial direction) 
we have determined front velocities, axial growth rates, and the wave numbers and frequencies of 
the unfolding vortex patterns with azimuthal wave numbers m = 0, ±1, respectively. 

PACS numbers: PACS nuniber(s): 47.20.-k, 47.54.+r, 47.32.-y, 47.10.+g 
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I. INTRODUCTION 

The Taylor- Couette system Q| of fluid flow in the annulus between concentric cylinders 
with the inner and the outer one rotating with different velocities is one of the simplest 
examples of a driven nonlinear dissipative system that shows spontaneous pattern formation 
out of an unstructured basic state that is stable at small driving 5]. This basic flow state 
is stationary and axially and azimuthally homogeneous and shows only a radial variation 
across the annular gap. It consists of a superposition of circular Couette flow (CCF) in 
azimuthal direction and of an annular Poiseuille flow (APF) in axial direction if as in our 



case an axial through-flow is imposed. Axially periodic vortex flow solutions bifurcate 
out of this homogeneous basic flow when the rotation rate of the inner cylinder is sufficiently 
high. These primary bifurcation thresholds to periodic vortex stuctures have been the aim 



of many linear stability analyses of the basic flow state 

For the radius ratio rj = 0.5 and the parameter ranges of rotation rates and through-flow 
investigated here in this work three spatiotemporally differing primary vortex structures are 
relevant: Rotationally symmetric, toroidally closed vortices with azimuthal wave number 
m = that move in downstream direction with the APF — for shortness we call this flow 
state Taylor vortex flow (TVF) although the presence of an axial through-flow modifles the 
genuine stationary TVF stucture. And, furthermore, spiral vortex flow (SPI) consisting of 
either left spiral vortices (L-SPI) with m = 1 or right spiral vortices (R-SPI) with m = —1. 

L-SPI and R-SPI are axial mirror images of each other in the absence of axial through- 
flow with the latter breaking the mirror symmetry of the former. While rotating azimuthally 
into the same direction as the inner cylinder L-SPI propagate axially opposite to R-SPI. 
This spiral dynamics is largely induced by the advective properties of the basic flow state. 
Furthermore, without through-flow the symmetry degenerate bifurcation treshold for these 
two symmetry degenerate SPI solutions is simultaneously also the bifurcation threshold for 
a vortex flow solution called ribbons 0]. This solution consists right at threshold of a linear 
superposition of L-SPI and R-SPI with equal amplitude and it becomes further away from 
threshold a genuine nonlinear vortex flow solution. However, here we are dealing only with 
linear vortex flow flelds that may be superimposed with arbitrary amplitudes as well as 
wave numbers and that are evolving separately from each other according to the linear fleld 
equations. Thus we do not need to discuss ribbons separately from our general investigation 
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of linear vortex modes with general axial and azimuthal wave numbers. 

In this work we quantitatively determine the influence of an axial through-flow on the 
spatiotemporal growth properties of linear perturbations of the basic flow state with az- 
imuthal wave numbers m = and m = ±1, i.e., of toroidally closed vortices and of spiral 
vortices, repectively. In each case we investigate (z) axially extended structures, (m) pulses 
of axially localized wave packets of vortices, and {in) vortex fronts. 

In Sec. |n]we describe the system, we briefly review the linearized Navier-Stokes equations 
(NSE) for the eigenvalue problem describing vortex perturbations, and we give details of our 
numerical procedure to solve the eigenvalue problem. In Sec. lIIII we discuss the spatiotempo- 
ral structure, symmetry properties, and bifurcation thresholds for onset of axially extended 
vortex perturbations of the form e'^i'^^+'^'p) with real axial wave number k and different az- 
imuthal wave numbers m in the absence and presence of an axial through-flow. In Sec. IIVI 
we consider axially localized wave packets consisting of superpositions of vortex eigenmodes 
of the linear NSE. Here we determine among others the boundary between convective and 
absolute instability of the basic flow against growth of vortices with a particular m by a 
saddle point analysis of the linear complex dispersion relation of the NSE over the plane of 
complex axial wave numbers. In addition we also determine the Ginzburg-Landau ampli- 
tude equation (GLE) approximation for the dispersion relation for the sake of comparison. 
In Sec. |V]we evaluate the spatiotemporal properties of linearly selected vortex fronts using 
a saddlepoint analysis of the dispersion relation. Also here we compare with GLE results. 
The flnal section contains a summary. 



II. SYSTEM 



Here we describe the system and we provide deflnitions and equations. Then we briefly 
review the linearized equations for the eigenvalue problem describing vortex perturbations 
of the basic flow state. Finally we give details of our numerical procedure to solve the 
eigenvalue problem. 
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A. Setup 

We consider the flow of an incompressible fluid in the annulus between two concentric 
cylinders of inner radius ri and outer radius r2 with a gap width d = r2 — ri. The boundary 
conditions at ri and r2 are no-slip. The angular velocity of the inner and outer cylinder is 
fli and respectively. The associated Reynolds numbers are 

Ri = -nili, R2 = -r^n^, (2.1) 
u V 

where v is the kinematic viscosity. An externally imposed axial through-flow is measured 
by the axial Reynolds number 

Re = ^{w) (2.2) 

where the mean axial velocity (w) averaged over the annular cross section describes the total 
through-flow. We use also the relative control parameters 

^~RURe) ' '~RiciRe = 0) ' 
measuring the relative distance of the inner Reynolds number Ri from the critical onset 
Ric of axially extended spiral vortices or Taylor vortices in the presence and in the absence 
{Re = 0) of through-flow, respectively In this notation 

is the critical threshold for onset of the vortex flow in question. The relation between /z and 
e is 

_ e - e,{Re) 

^-l + e^{Rey ^^-^^ 
With inflnitely long cylinders the only relevant parameter characterizing the geometry is the 
radius ratio 77 = ri/r2. 

The velocity fleld u of the fluid is described by the Navier-Stokes equations (NSE) for 
incompressible fluid flow 

9tu= V2u-i?i(u- V)u- Vp, V-u = 0. (2.6) 



Here and in the following we scale positions by the gap width d, the velocity u by the 
velocity riQi of the inner cylinder, time t by the momentum diffusion time d'^/u across the 
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gap, and the pressure p by prS^iu/d with p denoting the constant mass density of the fluid. 
Furthermore we decompose the velocity fleld 

U = UBr + Ve^ + WBz (2-7) 

into radial (m), azimuthal (f ), and axial (w) components using cylindrical coordinates r, </?, z. 
B. Basic flow state 

The basic flow state Uq that is reahzed in the absolutely stable regime of inner Reynolds 
numbers Ri below the thresholds for onset of Taylor and spiral vortex flow is rotationally 

symmetric, axially homogeneous, and constant in time. It consists of a linear superposition 
of circular Couette flow (CCF) in azimuthal direction, e^^, and of annular Poiseuille flow 
(APF) in axial direction, 

Uo = vccF{r)e^ + WAPFir)^^ (2.8) 
without any radial component. Here 

vccF{r)=Ar + B/r, (2.9) 



and 



with 



P + Cln{r) + D 
WAPF{r) ^ Re , (2.10) 



- n2/^i 



(2.11) 



{l-r])ln{7]) (l-r/)2' ^ ' ' 

E = (2.15) 
2 (1 -r7)2 1n(r7) ^ ' 
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C. Linear eigenvalue problem of vortex perturbations 

Let = {u^v^w^p) abbreviate the deviation fields from the basic flow state ()2.8|) . Then 
the general solution of the NSE linearized in the deviation fields can be written as a super- 
position of modes of the form 

V^(r, ^, z, t) = 0(r)e'('="+™'^)e"* (2.16) 

with axial wave number k = 27i/X and integer azimuthal wave number m. The complex 
amplitude functions 

0(r) = [U{r), V{r), W{r), P(r)] (2.17) 

depend on the mode indices k, m and the radial coordinate r. The characteristic exponent 
(t(/c, m) is in general complex. It is decomposed here as follows 

a = ^a + i'^a = -i -iuj (2.18) 

into the growth rate 7 and the characteristic frequency uo of the k — m mode. Substituting 
the above solution ansatz into the linearized NSE yields 



aU= { dl + -dr - ^^^ - k^ ]U + 2FV 



-drP - im { + FU ] - ikHU (2.19) 



aV ={d^^+ -a - -e\V + 2GU 



ly iy>^ 

im / 2 \ 

-—P + ini{—U-FV\-ikHV (2.20) 

aW =( dl + -9, - ^ - A:' W - %kP 
\ r r''' J 

-imFW -lU- ikHW (2.21) 

1 7771 

= d,U + -U + —V + ikW. (2.22) 

The solution of this eigenvalue problem yields the characteristic exponent a and the associ- 
ated eigenfunctions (f){r) as functions of k, m. Here 

F{r) = ^vccFir) , G = -R^A, (2.23) 
H{r) = RiWAPFir) , /(r) = drH{r) (2.24) 
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are quantities defining the basic flow state ()2.8p . The latter enters via the hnearized advective 
term of the NSE. 

In order to rewrite ()2.19I - I7.22|1 into a system of flrst-order differential equations - which 
is advantageous for numerical reasons - we introduce three additional complex amplitude 
functions llll 



X = drU + -U-P. 

r 

Y = drV + -V, 

r 

Z = drW. 



(2.25) 

(2.26) 
(2.27) 



Using ()2.25|) and the continuity equation ()2.22j) one can then eliminate the pressure in 
Eq. ()2.20j) by P = —X — ikW — ^V. All in all one obtains in this way a system of 6 coupled, 
first-order differential equations 

a^X = £ X, 



for the six variables 



with 



X={U,V,W,X,Y,Z)' 



(2.28) 



(2.29) 



C 



—i- 



-ik 
10 















1 



where 



L 2{i^-F) 
-2{i^ + G) L + ^ ^ -if 
J ^ L + e -ik 

r r 

2 

L = a + — + k'^ + imF + ikH. 



(2.30) 



(2.31) 
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D. Numerical procedures 

We have solved the eigenvalue equations ()2.28|) numerically with a standard shooting 
method subject to the six boundary conditions, 

U = V = W = at ri = r//(l - rj) and rs = 1/(1 - rj) , (2.32) 

which make the eigenvalue spectrum discrete. To integrate from ri to r2 we used a fourth- 
order Runge-Kutta formula with two step widths ( Ar = 1/200 and 1/400, respectively) 
for a Richardson extrapolation. A Newton-Raphson method [ld| was then used to find the 
roots of the complex determinant of the 3x3 matrix which ensures the vanishing of U, V, W 
at the outer cylinder. We therefore vary in the Newton-Raphson procedure only two of the 
parameters (a, r], i?2. Re, m, k) that enter into ()2.28|1 while keeping the others fixed. In 
this way we determine on the one hand the marginal threshold values of Ri and uj with 7 = 
for which the basic state is marginally stable against the growth of an extended perturbation 
with given m and real axial wave number k at specified parameters r/, R21 Re. On the other 
hand we calculate for given m, 77, i?2. Re the complex eigenvalue a over the complex wave 
number plane - including as special case also the real fc-axis. In each case we are interested 
only in the vortex modes with the largest growth rates for which the associated amplitude 
functions 0(r) display the least radial variation with the fewest number of nodes. 

We present here results for the radius ratio rj = 0.5 in a range of outer Reynolds numbers 
—150 < R2 ^ 50. In this parameter regime the vortex perturbations with the largest growth 
rates have in the absence of through-flow azimuthal wave numbers of either m = Oorm = ±l 



ll| . We investigate here linear properties of such vortices with m = and m = ±1 in a 
through-flow of Reynolds numbers < Re < 20. 



III. AXIALLY EXTENDED VORTEX STRUCTURES 

Here we discuss the spatiotemporal structure, symmetry properties, and bifurcation 
thresholds for onset of axially extended vortex perturbations with real axial wave num- 
ber k and different azimuthal wave numbers m in the absence and presence of an axial 
through-flow. 
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A. Spatiotemporal structure 

Structure and dynamics of the vortex modes ()2.16|) are dominated by the fact that their 
phases are constant on any cyhndrical surface, r = const, along hues given by the equation 

^o = -jV + ^Y^t. (3.1) 

Here the constant phase coming from the amphtude 0(r) is suppressed. Thus, on the (p — z 
plane of such an 'unrolled' cylindrical surface these lines of constant phase are straight. 



1. Taylor vortex like patterns — m=0 

For rotationally symmetric Taylor vortex like perturbations the line pattern of constant 
phases, kzo = n27c, is parallel to e,^. The m=0 pattern is stationary for Re = 0. Only for 
finite through-flow it propagates axially with phase velocity 

Wphase = ^ (3.2) 

that is proportional to Re. The main reason for this is that the azimuthal flow of the basic 
CCF state is precisely parallel to the vortex lines of constant phase while the APF flow 
being perpendicular to them can advect them. The latter happens in our axially periodic 
system that does not exert any phase pinning at the axial boundaries as soon as Re > 0. 



2. Spiral patterns — m / 

The vortex modes ()2.1(i|l with axial wave number m ^ have spiral structure. When 
m/k is positive (negative) the lines of constant phase ZQ{(f, t) (jH.lll wind in a left spiral L-SPI 
(right spiral R-SPI) around the cylindrical surface r = const with negative (positive) slope 
d,pZQ = —m/k. The lines of constant phase and with it the whole spiral stucture rotates in 
if rigidly with angular velocity 

= — • (3.3) 

m 

In the absence of an externally imposed through-flow, Re=0, this rotation proceeds for L- 
SPI and R-SPI alike into the same direction as the rotation of the inner cylinder. The reason 
is that the spiral perturbations are advected by the inner part of the azimuthal CCF which 
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is relevant for the centrifugal instability leading to vortex generation. A model explaining 
this effect is presented in |l7|. 

There are two immediate consequences of this advective origin of the spiral dynamics 
induced by the inner cylinder's rotation: (i) With the latter being by definition positive - 
in this work the inner cylinder is taken to rotate in positive yj-direction - also uj{k, m)/m is 
positive for Re=0. Hence, say, an m = 1 (m = — 1) spiral has positive (negative) frequency 
for Re=0. (ii) Then a L-SPI (R-SPI) being defined hj m/k > {m/k < 0) propagates for 
Re=0 upwards (downwards) with positive (negative) axial phase velocity 

Wphase = ^ = '^'^^P^ ^^-^f 

that is directly related to its positive angular velocity ipspi- 

An externally applied axial through-flow changes the axial phase velocities and frequencies 
of the m = and m ^ vortex modes roughly proportional to Re, i.e., 

Wphase{k, m, Re) - Wphase{k, m, Re = 0) (x Re. (3.5) 

Simultaneously the rotation rates, 0sp/ = Wphasek/m, of the spirals are changed accordingly. 
Thus, in each case the vortex frequencies are largely determined by the basic state's advective 
properties, i.e., by the combination of azimuthal advection by vccf and axial advection by 
wapf- 



B. Symmetries 

Here we consider symmetry properties of axially extended vortex perturbations of the 
form ()2.16|) with real wave number k. Symmetry relations between different vortex fronts 
with complex wave number Q are discussed later on in Sec. IVB 21 

Table IJ shows the symmetry transformations that leave the eigenvalue problem un- 
changed. They reflect (i) that C transforms under complex conjugation (indicated by an 
overline) as 

C{k,m,a) = C{—k, —m,a) (3.6) 

and (a) that the NSE ()2.6|) are invariant under an axial reflection {z,Re,w) — > 
{—z, —Re, —w). 
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Thus one infers, for example, that the growth rate (frequency) of the characteristic expo- 
nent a = 'J — iuj for m = vortices is an even (odd) function of k and Re. For perturbations 
with m 7^ one finds that 

'j{k,m, Re) = 'y{~k, —m, Re) = 'y{—k,m,—Re)='y{k,—m,—Re) (3.7) 
u{k,m, Re) = —uj{—k, —m, Re) = uj{—k,m,—Re) = —uj{k,—m,—Re) (3.8) 

and that the spatiotemporal structure including the amplitude functions of a L-spiral per- 
turbation {k/m > 0) at Re > is the same as that of a R-spiral {k/m < 0) at Re < 0. 
Note, however, that any finite through-flow breaks the axial mirror symmetry between L- 
and R-spirals at Re = so that, among others, 

-fni-Re) = ^L{Re) ^ iR{Re) = -fL{-Re) (3.9) 

when Re ^ 0. But the symmetry relations are such that it suffices to investigate, say, 
positive k combined together with either (i) m > only for positive and negative Re or, 
equivalently, {ii) m positive and negative for Re > only in order to get the complete linear 
information on both spiral vortex types. 

C. Bifurcation thresholds 

Figure n] shows the critical bifurcation thresholds Ric{R2) for rn = and m = ±1 vortex 
patterns with the respective critical wave numbers, kc{R2), as functions of R2 in the absence 
of through-flow. 

The vertical lines in Fig. Q mark the two outer Reynolds numbers R2 = and R2 = 
—125 for which we show in Fig. El as representative examples how the critical thresholds 
evolve with through-flow Reynolds number Re. The above discussed symmetry relation 
7ij(— -Re) = 'jL^Re) between the growth rates of R- and L-spirals implies the corresponding 
relation between the respective bifurcation thresholds (full and dashed lines in Fig. |21). In 
the remainder of this paper we therefore restrict ourselves without loss of information to 
positive Re. 

For small Re > the axial flow stabilizes the basic state against growth of TVF (m = 0) 
and R-SPI (m = —1) perturbations. On the other hand, the bifurcation threshold for L- 
SPI (m = 1) vortex patterns that propagate into the same direction as the through-flow 
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decreases at small Re and increases only at larger Re. The upwards shift of the m = 
threshold with increasing Re is stronger than the one for m = 1. Thus, eventually the latter 
comes to lie below the former and consequently the growth of L-spirals propagating into the 
same direction as the through-flow is favored for sufficiently large Re even when R2 = 0. 

Dotted lines in Fig. IHl show the reduced critical threshold curves ec(-Re) ()2.4|1 as functions 
of Re for the two representative outer Reynolds numbers R2 = and R2 = —125. The other 
lines in Fig. El are discussed in Sec. IIVI Our numerical results for e^Re) that were obtained 
in steps of 6Re = 1 were fitted in the range Re = —20 ■ ■ • 20 to the following expression 

f = aiRe + a2Re^ + a^Re^ + a^Re\ (3.10) 

The fit parameters a„ are listed for different R2 in Tables |n] and IIIII for TVF (m = 0) 
and L-SPI (m = 1), respectively. The threshold curves for R-SPI (m = —1) are obtained 
according to Sec. IIII Bl from those for L-SPI by Re —>■ —Re, i.e., by changing the sign of the 
odd coefficients in Table UTTl 

IV. LOCALIZED VORTEX PERTURBATIONS 

So far we have considered axially extended vortex perturbations of the form iplr, z, t) = 
^(^^^^i{kz+m<fi) ^\^[q[^ are single eigenmodes of the operator C ()2.30p . For supercritical control 
parameters, -Ri > Ric{m), a finite band of axial wave numbers can grow and with it also 
axially localized wave packets consisting of superpositions of vortex modes. 

A. Vortex packets 

Let us consider first an infinitesimal initial perturbation with azimuthal wave number m 
that is axially localized, i.e., a wave packet that consists of a superposition of vortex modes 
of different k but common m - an initial perturbation containing different m-modes would 
be just a sum of the above described vortex packets that would evolve independently of each 
other as long as the linear description is valid. 

After fast transients have decayed a pulse like perturbation survives with axial wave 
numbers within the unstable band centered around the wave number of maximal growth 
kmaxiiT^) — kc{m). Siucc the above described wave packet contains vortex modes that can 
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grow the pulse will grow as well. Simultaneously the pulse center travels axially for small e 
with the critical goup velocity 

.„=^ . (4.1) 



-"9 



dk 



Hence when at a bifurcation threshold the frequency is nonzero with a finite group velocity 
then the supercritical spatiotemporal growth behavior of an axially localized vortex pertur- 
bation differs significantly from an axially extended vortex mode. The growth of the latter 
is axially uniform which is not the case for the former. 

Furthermore, one has to distinguish between two different supercritical regimes for the 
former: (z) In the so-called convectively unstable parameter regime of the basic state the 
vortex packet moves with the velocity Vg faster away than it grows - while growing in 
the frame comoving with Vg the pulse moves out of the system so that the basic state is 
restored jisl . llQ^ . In other words, the two fronts that join the pulse intensity envelope to 
the structureless state propagate both in the direction in which the packet center moves. 
(ii) In the so called absolutely unstable parameter regime the growth rate of the packet is 
so large that one front propagates in the laboratory frame opposite to the center motion. 
ThuSj_the packet expands not only into the direction of the pulse motion but also opposite 
to it [18] so that eventually the initial perturbation can fill the entire system. But the linear 
growth analysis of the vortex fields does not determine in what nonlinear final stable state 
the system will end nor what possible intermediate nonlinear transient behavior might occur. 

However, this analysis has an important implication for experiments with through-flow: 
Developed vortex structures can be seen in finite systems with vortex suppressing inlet 
conditions only in the absolutely unstable regime which is typically realized at larger Ri > 



Ric (cf. Fig. El) if one leaves aside noise- sustained patterns [20|, |2l| in the convectively 
unstable regime. In this latter regime the vortex front that connects to the zero amplitude 
inlet condition moves downstream (we assume that the fronts of our forwards bifurcating 
vortex structures are linearly selected thus excluding the buildup of nonlinear fronts that 
might revert their propagation direction). In the absolutely unstable regime, on the other 
hand, an upstream motion of this front is stopped by the inlet condition at a characteristic 
downstream growth length from the inlet. This growth length of the downstream evolving 
vortex structure diverges when approaching the boundary between convective and 

absolute instability from the latter regime. 
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It is remarkable that based on the experimental observations of Takeuchi and Jankowski 
0] such a behavior was discussed already in 1979 [cf. figures 4 and 5(a) and their discussion 
in Sec. 6 of ref. j^], albeit without invoking the concept of absolute and convective instability 
1^1 which was introduced to a broader fiuid dynamics community only a few years later 

B. Boundary between convective and absolute instability: saddle point analysis 

The boundary between convective and absolute instability of the basic fiow against growth 
of vortices with a particular azimuthal wave number m is marked by those parameter combi- 
nations for which one of the fronts of the linear packet of m- vortices reverts its propagation 
direction in the laboratory frame: In the convectively unstable parameter regime this front 
propagates in the same direction as the center of the packet, in the absolutely unstable 
regime it moves opposite to it, and right on the boundary between the two regimes the front 
is stationary in the laboratory frame. 

This parameter combination can be determined by a saddle point analysis of the linear 
complex dispersion relation a{Q) over the plane of complex wave numbers [l^ 

Q = ^Q + i'^Q = k-iK . (4.2) 

Here we do not display the dependence of 

(t(Q) = 3fta(g) + iQaiQ) = 7(Q) - zcu(Q) (4.3) 

on the parameters Ri,R2, Re, t] and we also do not indicate that the dispersion relations for 
vortex perturbations with different azimuthal wave numbers m are different. The boundary 
condition of no growth for a front that is stationary in the laboratory system is 

3fJa = (4.4) 

at the appropriate saddle position, Q*, of a{Q) in the complex wave number plane Here 
Q* follows from 

- 0. (4.5) 



dQ 

For fixed m, R2, Re,ri eqs. (l4.4l4.5|) yield Ric-a- Here and in the following the subscript 
c — a identifies boundaries between convective and absolute instability. Thus, e.g., the basic 



February 2, 2008 



15 



flow state is convectively [absolutely] unstable against vortex perturbations with azimuthal 
wave number m for Ridm) < Ri < Ric-aim) [Ri > Ric-a{^)] and absolutely stable when 
Ri < Ric{m). 

We are interested here in the Re dependence of these thresholds and we will discuss to 
that end the reduced boundary quantities 

^-^"^'^ = -RUW - ' ^'-'^ 

and 

-Ri - (Re) 

eJRe) = — — r — 1, where x = coix = c — a. (4.7) 

Ric[Re = 0) 

Because of the symmetries of cr{Q), of its saddle-point Q*, and of the resulting boundary 
Ric-a{rn, Re) one has 

(4.8) 
(4.9) 
(4.10) 
(4.11) 

so that it suffices again to investigate only positive through-flow. 



k*{—m, 


-Re) 


= k*{m,Re) 


K*{-m, 


-Re) 


= -K*{m,Re) 




-Re) 


= fic-a{m,Re) 




-Re) 


= ec-a{m,Re) 



C. Numerical procedures 

In order to determine the boundary between convective and absolute instability via the 
solution of eqs. ()4.4l 14. 5|) one has to evaluate the dispersion relation cr(Q) for complex Q 
To that end we solved the eigenvalue problem ()2.28j) of the full field equations for 
complex Q (cf., Sec. lIV(T2|l . In addition and for comparison we used for a{Q) the Ginzburg- 
Landau amplitude equation approximation which only requires knowledge of o"(A;) along the 
real fc-axis. 
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1. Ginzburg- Landau amplitude equation approximation 

Within this approximation the dispersion relation a{Q) for vortex modes is expanded in 
Q and -Ri around the critical point kc, Ric 

+ (^) (^^^^^c) + h.o.t. (4.12) 

= -iujc - iVgiQ - kc) ^(1 + ici){Q - kc)"^ 

+ ^i±^/. + /..oi.. (4.13) 

To 

The expansion coefficients in the above expressions appear also in the linear parts of the 
Ginzburg-Landau amplitude equation P]. They are obtained from the numerical solution of 
the eigenvalue problem ()2.28|1 of the full field equations for real k close to the critical point 

Within the GLE approximation one gets from eqs. ()4.4l 14. 5|) 

k* K* 
2 2 



and 



e,_,(i?e) = ec{Re) + ' \ [1 + ee(i?e)] . (4.16) 

Note that the GLE coefficients fc,, tq, Cq, Ci depend on m, i?2, -Re, rj. For m = and -R2 = 
also the nonlinear coefficients of the GLE have been obtained for several 7] as functions of 
Re \m. 



2. Dispersion relation a{Q) of the NSE for complex Q 

To assess the quality of the GLE results for the boundary between convective and absolute 
instability we determined the dispersion relation of the NSE not only for real k close to the 
critical point but for complex axial wave numbers Q that lie in the vicinity of the relevant 
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saddle locations Q* of cr(Q). Having determined <7{Q) as described in Sec. IIIDI we solved 
the equations ()4.4|4.5|) for Ric-a in the form 



0; 



OK 



0; 



dl{Q) 



dk 



(4.17) 



that follows from using the Cauchy-Riemann relations for cr{Q). 



D. Results 

In Fig. m we show the i?e- variation of the real and imaginary parts of the relevant saddle- 
point Q* = k* — iK* at the boundary Ric-a for two characteristic cases: (a) L-SPl (m = 1) 
at -R2 = and (b) R-SPI (m = — 1) at i?2 = ^125. In each case full (dashed) lines were 
obtained from the correct NSE (approximate GLE) dispersion relation cr{Q). 

Case (a) is representative for a situation where the GLE approximation reasonably well 
reproduces the correct result - at least for small Re - and starts to deviate significantly 
only for larger Re. On the other hand, in case (b) the GLE approximation to k* displays 
a smooth variation with Re that reflects the smooth variation of the saddle Qgle ^4.14|1 
of o-gle{Q) f!4.13|) while the real part k* of the saddle location of the correct dispersion 
relation undergoes a dramatic change around Re ~ 0.8. The reason is that the surfaces 
of 7((5) = 3?cr(Q) over the Q-plane for the two largest eigenvalues intersect and change 
their order along the 7 axis. Thus, the saddle that is relevant for the c — a transition and 
that has the largest 7 value switches from one eigenvalue surface to another. Similarly the 
eigenvalue surface might have another saddle and their ordering in 7 changes. In contrast 
to that the GLE approximation produces only a single eigenvalue tracing out the smooth 
surface (Jgle{Q) ()4.13|) . 

To give an impression of such an intersection of the correct dispersion surfaces we show 
in Fig. El their real parts over the complex Q— plane at i?e = 5, = Hc-a- Full lines 
in Fig. ini show different constant-fc sections through them in the intersection range. At 
Re = 5, fi = fic-a the saddle has moved already away from the intersection region. The 
saddle coordinates {Q* ^ 4.73 — il.76,7* = 0) are indicated in Figs. IH El by full dots. 

While the saddle locations Q*gle ^^e GLE approximation acLEiQ) ()4.13p can differ 
substantially from the saddle Q* of the correct dispersion relation the difference in the 
boundaries between convective and absolute instability is typically much less pronounced 
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- cf, e.g., the dashed and full curves for ec-a in Fig. El There the GLE results (dashed 
lines) agree in each case quite well with the correct boundaries (full lines) up to, say Re = 
5. However, as a consequence of the typical increase of the reduced boundary /ic-a with 
increasing Re the quality of the GLE results for ^c~a generally deteriorates. After all the 
GLE is strictly valid only for — > 0. Fig. [3 shows an example (m = 1, i?2 = —125) where 
the GLE prediction for the boundary shows even a qualitative different variation with Re 
for, say. Re > 15. We have observed similar behavior - partly at larger Re - also for other 
combinations of m, R2. 

Tables HTl and ITTTl contain the coefficients for fitting our results for ec-a(^ = 0,R2,Re) 
and ec-a(^ = 1,-^2, -Re) in the same way as described in Sec. IIII (^1 Also here the boundary 
curves for R-SPI (m = —1) are obtained according to Sec. IIII Bl from those for L-SPI (m = 1) 
by Re —Re, i.e., by changing the sign of the odd coefficients in Table ITTTl 



V. FRONTS AND PULSES 



The vortex fronts that we are investigating here and that appear also as constituents of 
vortex pulses are perturbations of the basic state where the fields (locally) have the form 

in the laboratory frame. This form describes long-time spatiotemporal properties of uniquely 
selected linear fronts [5]. Axial front velocity f *, axial growth rate K* of the front envelope, 
wave number k* of the vortex pattern under the front, and its frequency uj* are determined by 
the saddle behavior of the linear complex dispersion relation cr((5, to) = '~f{Q, m) — iuj{Q, m) 
of the field equations over the plane of complex axial wave numbers Q = k — iK . Since 
superpositions of vortex perturbations with different m evolve independently from each other 
within the linear description as mentioned already in Sec. IIV Al we consider here only fronts 
of vortex perturbations that have a common azimuthal wave number m and we do not always 
display the latter explicitly. 
The saddle condition is 

And the stationarity requirement that the temporal growth rate of the front vanishes in the 
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frame comoving with the front velocity v* demands that 



= gfj [(r{Q) + ivQ] 



-f{Q*) + v*K*. 



(5.3) 



Q=Q* 



We combine fj5.2|) and ()5.3|) into the three equations 



* 

V = 



K* 



dK Q* 



dk Q* 







(5.4) 



that we have solved for w*, A;*, K* . 

A. Notation 

The front envelope of ()5.H) varies axially with e^*^. If K* > 0, then the perturbation 
grows at z = — oo out of the basic state. We call such a front to be of type + and identify the 
associated front properties by a subscript +. On the other hand, for K* < we have a front 
of type — with an intensity envelope that joins at 2; = oo with the basic state. So the two 
subscripts ± identify the axial variations of the front envelopes. A pulse-like perturbation 
of the basic state would consist suffiently away from its center of a + front for 2; ^ — cxd and 
of a — front for z ^ 00. Correspondingly the eqs. ()5.4|) have two different solutions: one 
of them describes a — front and the other one a + front. A schematic plot of the different 
envelope types can be seen in Fig. |H1 

We identify the vortex pattern that is unfolded under a front either by its azimuthal wave 
number m or by the superscripts T, L, R. Here T refers to a TVF-like pattern of vortices 
with m = 0, L denotes L-spiral vortices with m = 1, and R identifies R-spiral vortices with 
m = —1. Hence m or equivalently the superscripts identify the spatiotemporal structure of 
the vortex pattern growing under the front. In this work we restrict ourselves to these three 
vortex varieties. Since they appear under the above described two front envelopes there are 
six different fronts that we have investigated here. 

B. Results 

In Fig. ini the results of our investigations are shown for fronts with azimuthal wave 
numbers m = 0, 1, and -1 for axial through-flow Reynolds numbers Re =0, 10, and 20. 
In each case the outer cylinder is at rest, R2 = 0. The front properties are presented as 
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functions of fi. Within each 2x2 block of figures (a)-(i) the left column shows real part k* 
and negative imaginary part K* of the saddle point and the right column shows the front 
velocity v* and the frequency u*, respectively. They all start at /i = since vortex growth is 
possible only above the critical threshold, i.e., for Ri > Ric{Re). The two critical fronts at 
this threshold are degenerate (fcl = k*^ = kc, uj'L = uj'^ = uJc, f * = = Vg) with vanishing 
axial growth rates = = 0. For /i > 0, however, + and — fronts differ from each 
other thus reflecting differences in the respective saddle points. 

The variation of k*, K*, v*, and u* with /i and Re is best understood by comparison with 
the corresponding GLE approximation (cf. Sec. IV B Hi and by invoking symmetry relations 
(cf. Sec. IVB 2|1 . We therefore continue the discussion of our results in the following two 
sections IYRTI and IVR2l 



1. Ginzhurg- Landau amplitude equation approximation 

The saddle point analysis ()5.2I - IK^ of the GLE approximation ()4.1Hj) yield the following 
front properties 



kl = k,-c^Kl, (5.6) 
vl = t;,-2(l + c?)^K;, (5.7) 

To 

7; = -v*^K*^, (5.8) 

= Uc + Vg{k^ - kc) + (ci - cq) — . (5.9) 

To 

for the two fronts with K*^ > and i^* < 0, respectively. All quantities appearing in (j5.5l - 
I5.9|l depend on whether they refer to T (m = 0), L (m = 1), or R (m = -1) vortex fronts. 
These GLE results are shown in Fig. El by dashed lines. 

They reasonably well describe the small-/i behavior of the correct front properties (full 
lines in Fig. E)) which were obtained from the correct dispersion relation of the NSE: As 
predicted by the small-/i GLE approximation (j5.5l - l^^ one finds that for small fi the axial 
growth rates K* vary oc ^/I) that consequently also k* — kc and v* — Vg vary oc y/JI, and 
that u* — UJc can have in addition also a contribution oc /i when ci — cq 7^ 0. The latter is 
the case for m = ±1 irrespective of Re and for m = if Re 7^ 0. 
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Note that for TVF fronts with m = 0, Re = the GLE predicts oo* = whereas the correct 
dispersion in Fig. Efa) seems to show for small fi a variation of uj* oc that is beyond the 
range of applicability of the GLE. Thus, under the linear part of a moving TVF front there 
should be a non zero phase propagation in the laboratory frame with phase velocity uj*/k*. 
The analogous behavior was found also for convection rolls in the Rayleigh-Benard system 

El- 

At the boundary between convective and absolute instability, fi = fj^c-a, the velocity v*^ 
and the temporal growth rate 7^ of the + front vanish in the laboratory frame while is 
positive there (and given by 2 Vg within the GLE approximation). Note that according to 

Fig. Eland tables ITTland lllll /i^ is zero for m = 0, Re = and very small for m = ±1, Re = 0. 

Only for large enough through-flow fic-a becomes sizeable. In the convectively unstable 
regime < ^ < ^c~a both, the + front as well as the — front of a vortex pulse move into the 
same downstream direction as the through- flow, Q < v*^ < v*_. In the absolutely unstable 
regime /i > /Xc-a, however, the + front moves upstream and the —front moves downstream, 
v\ <Q <v* . 

In ref. 0| it was remarked that the phase velocity of spiral patterns in axial flow through 

deviates from the critical one, cUc/fcc? of axially extended 



a system of radius ratio r] = 0.95 
vortex perturbations. While we have done only calculations for rj = 0.5 our results of Fig. |H1 
and of the small-/! GLE approximation ()5.5I - shed some light on the existence of such 
deviations. Since the experimental vortex structures grow in downstream direction under 
intensity fronts one strictly speaking would have to compare with the phase velocity under 
such fronts that connect in the absolutely unstable regime to the fully developed downstream 
vortex pattern. Such a nonlinear analysis has been done for m = patterns |23i] but not for 
spirals. However, already the phase velocities uj*^/k*^ under our linear + fronts that grow in 
downstream direction differ from the corresponding critical phase velocities Uc/ of axially 
extended patterns — cf. Fig. 121 and eqs. ()5.5l - l5^ . 

2. Symmetries 

The front properties shown in Fig. IHl are largely influenced by the symmetry properties 
of the system without through-flow although a finite Re changes them. 

Invariance of the field equations under z —>■ —z for Re = implies that stationary 
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perturbations with Uc = (TVF) under a + front are mirror images of those under a — 
front. This imphes for Re = the symmetry relations 

{K, uj, t;)f = -{K, 00, vyj- fcf = k*^. (5.10) 

These symmetry properties of the two front types of T-perturbations can be seen in Fig. 
Ha). 

Now consider L and R perturbations. Here the invariance of the field equations under 
z —>■ —z for Re = implies first of all that a spatially extended L-SPl with uniform amplitude 
is the mirror image of a spatially extended R-SPI. Furthermore, a L-SPI under a + front 
with positive K is symmetry degenerate with a R-SPI under a — front with negative K. 
Similarly a R-SPI under a + type front is the mirror image of a L-SPI under a — front. This 
implies for Re = the symmetry relations 

{K,uj,v)l'' = -{K,uj,vyJ^; kf = k*_'^, (5.11) 



{K, uj, v)f = -{K, u, vy_^; kf = k*^. (5.12) 

They can be seen to be realized in Fig. IHl^b) and Fig. Ofc). 

The GLE approximation ()5.5|) - ()5.9|) shows beyond the relations ()5.10p - ()5.12|) the follow- 
ing additional relations 

kf = k*^ = kj ; u*^ = 00*^ = (5.13) 
j^*L ^ _^*L ^ j^*R ^ _K*_R^ (5.14) 

ujf = u*^ = -uf = -uj*^ (5.15) 

that follow from the fact, that = = Cq = cf = 0, k^ = k^, = —Vg, = Tq, 
= -c^, cf = -cf , and = for i?e = i26|J. 

VI. SUMMARY 



We have determined the influence of an axial through-flow on the spatiotemporal growth 
behavior of structurally different vortex perturbations of the basic Couette-Poiseuille flow 
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in the Taylor- Couette system with radius ratio rj — 0.5. To that end we have solved the 
linearized NSE numerically with a shooting method for vortex perturbations with azmuthal 
wave numbers m = (TVF), m = 1 (L-SPI), and m = — 1 (R-SPI) in a wide range of the 
parameters Re,Ri, and i?2- Here symmetry properties allowed us to restrict ourselves to 
positive through-flow Reynolds numbers Re. For each of the three different vortex varieties 
we have investigated {i) axially extended vortex structures with homogeneous amplitudes, 
(ii) axially localized vortex pulses consisting of a linear superposition of axially extended 
vortex modes with different real axial wave numbers k, and {in) vortex fronts. 

Central to our analysis is the determination of the complex dispersion relations a{Q) of 
the linearized NSE for vortex modes with the three different m. We have evaluated a over 
the plane of complex wave numbers Q = k — iK for patterns (iz, in) and along the real 
/c-axis for pattern (i). We have also determined the Ginzburg-Landau amplitude equation 
approximation acLEiQ) in order to analyze its predictions for the vortex stuctures {ii, Hi) in 
comparison with the correct NSE dispersion relation a{Q). In each case symmetry relations 
are elucidated. 

First we have evaluated the critical bifurcation thresholds Ric{Re, R2,m) for axially ex- 
tended vortex structures. Then using a saddle-point analysis of cr[Q) we have determined 
the boundaries Ric-a{Re, R2,'m') between absolute and convective instability of the basic 
state at which one of the fronts of the expanding vortex pulses reverts its propagation direc- 
tion in the laboratory frame. Here we have elucidated also in some detail how the different 
saddle topologies of cr{Q) and of aaLsiQ) explain some of the shortcomings of the latter. 
Fit parameters for power-law expansions of the reduced boundaries Cc, ec-a, and e^^^ up to 
Re"^ are listed in two tables. 

Finally we have determined the linearly selected front behavior of growing vortex patterns 
with m = 0, ±1 for i?2 = under two different types of front intensity envelopes: type -|- 
shows growth in positive z-direction while type — locates growth in negative z-direction. 
The combination of the three different dynamics of the constituent vortex modes (m = 
0, ±1) and of the two different spatial intensity profiles (+, — ) leads to six different fronts 
^ ^t{k*z-uj*t)^K*{z-v*t) laboratory frame. Their velocity v* , spatial growth rate K* , 

wave number A;*, and frequency uj* as determined via a saddle point analysis of the respective 
dispersion relations differ in general from each other in the presence of a through-flow. 
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FIG. 1: Critical bifurcation thresholds i?ic(i?2) for m = and m = ±1 vortex patterns with the 
respective critical wave numbers, kc{m, R2), as functions of R2 in the absence of axial through-flow. 
The vertical lines mark the two representative outer Reynolds numbers R2 = and R2 = —125 
that are investigated in more detail in this work. The radius ratio is rj = 0.5. 
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FIG. 2: Evolution of the critical bifurcation thresholds Ric{Re) for m = and m = ±1 vortex 
patterns with through-flow Reynolds number Re. The two outer Reynolds numbers R2 = and 
R2 = —125 are marked in Fig. ^ The radius ratio is = 0.5. 
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FIG. 3: Stability boundaries of the basic flow state as functions of the through- flow Reynolds 
number Re. Dotted lines show the reduced critical bifurcation thresholds Cc (|2.4() for axially 
extended vortex patterns. The full (dashed) boundary lines, ec-a, between the convectively and 
absolutely unstable parameter regions for vortex growth were obtained from the eigenvalues of the 
full NSE (the GLE approximation) - cf. Sec. IIVI The radius ratio is r/ = 0.5. 
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FIG. 4: Evolution of the coordinates k* and K* of the saddle point that determines the boundary, 
Ric~a, between convective and absolute instability with through-flow Reynolds number Re. Full 
(dotted) lines are evaluated with the correct NSE (approximate GLE) dispersion relation. Filled 
circles in (b) at Re = 5 mark the coordinates of the saddle in Fig. \^ The radius ratio is tj = 0.5. 
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FIG. 5: Real parts 7(Q) of the two biggest eigenvalues of the NSE over the complex Q— plane for 
m = —1, i?2 = —125, Re = 5, ^ = jic-a, V = 0.5. The filled circle marks the saddle at ^{Q*) = 
that determines the boundary beteen convective and absolute instability. 
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FIG. 6: Full lines are sections through the surfaces of ^{Q) (Fig. IS} for the two biggest eigenvalues 
of the NSE at constant values of k in the vicinity of the intersection of the surfaces. Dashed 
lines show j{K) obtained from the GLE approximation to the dispersion relation. Parameters are 
m = —1, R2 = —125, Re = 5, /i = /ic-ai 1] = 0.5. 
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FIG. 7: Stability boundaries of the basic flow state as functions of the through-flow Reynolds 
number Re for m = 1, i?2 = —125. The full (dashed) boundary lines, ec-a, between the convectively 
and absolutely unstable parameter regions for vortex growth were obtained from the eigenvalues 
of the full NSE (the GLE approximation) - cf. Sec. IIVI Dotted lines show the reduced critical 
bifurcation thresholds Cc (|2.4p for axially extended vortex patterns. The radius ratio is 77 = 0.5. 
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FIG. 8: Schematic plot of different vortex fronts. Thick full lines in (a) and thin full lines in (b) 
show intensity envelopes of + type and of — type fronts connecting to the basic state at z ^ —oo 
and at z — > oo, respectively. This line convention — thick ones for + type fronts and thin ones for 
— type fronts — is used also in Fig. |^ Dashed lines indicate the vortex field growing under the 
front. 
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FIG. 9: Front properties of vortices with azimuthal wave numbers m = 0, 1, —1 in systems with 
and without through-flow Re for = 0, = 0.5 as functions of /x = R\l R\c{Re) — 1. Within each 
2x2 block of figures (a)-(i) the left column shows the axial wave number A;* and the axial growth 
rate K* . The right column shows the front velocity v* and the frequency in the laboratory 
frame. Thick lines correspond to + fronts, thin lines to — fronts, respectively. Full (dashed) lines 
result from the saddle point analysis of the dispersion relation of the NSE (GLE approximation). 
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Re 
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Re 
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U ^ 
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w - 


-w 


-w 


X 


X 


X 


X 


X 


Y 


Y 


Y 


Y 


Y 


z 


Z 


Z - 


-Z 





TABLE I: Transformation behavior of the eigenvalues and eigenfunctions of the eigenvalue problem 
1)2.281 - 12.31(1 under symmetry operations. Here A denotes the identity, B complex conjugation, 
C axial reflection, and D complex conjugation (indicated by an overbar) combined with axial 
reflection. 
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R2 


-150 


-125 


-100 


-75 


-50 


-25 





25 


50 




a2 * 10^ 


0.590 


0.795 


1.089 


1.502 


2.083 


3.115 


3.679 


2.447 


1.307 


a4 * 10^ 


5.181 


3.966 


1.151 


-2.804 


-3.385 


-15.62 


-34.97 


-14.55 


-0.958 




a2 * 10^ 


1.181 


1.451 


1.854 


2.516 


3.906 


5.588 


6.238 


4.338 


2.483 


a4 * 10^ 


0.116 


0.139 


0.182 


0.293 


-2.583 


-4.642 


-5.716 


-3.682 


-1.740 




a2 * 10^ 


1.161 


1.416 


1.800 


2.517 


4.068 


6.547 


7.732 


5.152 


2.762 


04 * 10^ 


-0.164 


-0.268 


-0.481 


-1.011 


-2.106 


-3.807 


-4.900 


-3.022 


-1.425 



TABLE II: Fitparameters for the -Re-dependence of the stabihty boundaries of the basic state 
against growth of TVF (m = 0) perturbations for different R2: reduced critical bifurcation thresh- 
old ec of axially extended vortex patterns and boundary ec-a between convectively and absolutely 
unstable parameter regime. Here e^^ is obtained from the GLE approximation (cf. Sec. IIV C 1|) . 
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R2 


-150 


-125 


-100 


-75 


-50 


-25 





25 


50 


Ec = aiRe + a2Re^ + a^Re'^ + a^Re'^ 


ai * 10^ 


-1.582 


-2.365 


-3.513 


-4.887 


-6.143 


-5.833 


-3.372 


-1.402 


-0.514 


02 * 10^ 


0.957 


1.199 


1.496 


1.925 


2.520 


3.174 


3.197 


1.806 


0.581 


as * lO'^ 


-9.181 


-9.001 


-8.485 


-7.678 


-6.741 


-4.172 


-3.832 


-3.572 


-3.685 


a4 * 10^ 


1.018 


-0.378 


-0.180 


-0.563 


-3.673 


-21.77 


-36.15 


-5.707 


19.92 


ec-a = oo + aiRe + a2Re^ + a-^Re"^ + a^Re'^ 


ao * 10^ 


2.358 


5.497 


3.775 


-1.237 


4.095 


12.25 


18.96 


10.09 


2.610 


ai * 10^ 


-3.268 


-3.578 


-3.048 


-0.449 


0.926 


1.159 


1.012 


1.217 


1.455 


02 * 10^ 


1.532 


1.964 


2.618 


3.539 


4.351 


5.535 


6.036 


4.357 


2.603 


as * 10^ 


8.408 


11.87 


14.97 


9.992 


6.487 


4.919 


3.680 


1.367 


-0.143 


a4 * 10^ 


-0.169 


-0.368 


-0.815 


-1.977 


-2.809 


-4.068 


-4.827 


-3.203 


-1.549 


e'^^^ = qq + aiRe + a2Re^ + osi^e^ + a^Re^ 


ao * 10^ 


1.101 


-2.541 


1.405 


2.529 


3.458 


3.998 


3.641 


1.320 


0.062 


ai * 10^ 


-4.117 


0.755 


-1.044 


-0.318 


0.831 


2.063 


1.495 


1.377 


1.625 


a2 * 10^ 


2.030 


2.792 


2.820 


3.435 


4.437 


6.166 


7.265 


4.981 


2.723 


as * 10^ 


-8.845 


-30.34 


-7.121 


4.780 


10.58 


8.229 


-8.191 


-12.58 


-10.03 


a4 * 10^ 


-0.822 


-2.713 


-1.136 


-1.004 


-1.674 


-2.923 


-3.799 


-2.234 


-0.854 



TABLE III: Fitparameters for the i?e-dependence of the stabihty boundaries of the basic state 
against growth of L-SPI (m = 1) perturbations for different R2: reduced critical bifurcation thresh- 
old ec of axially extended vortex patterns and boundary ec~a between convectively and absolutely 
unstable parameter regime. Here is obtained from the GLE approximation (cf. Sec. IIV C 1() . 
The results for R-SPI (m = — 1) perturbations are obtained according to Sec. IIIIBl from those for 



